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ABSTRACT 



We report spectra obtained with the Spitzer Space Telescope in the A = 14 
- 35 /im range of 19 nearby main- sequence stars with infrared excesses. The six 
stars with strong dust emission show no recognizable spectral features, suggest- 
ing that the bulk of the emitting particles have diameters larger than 10 /j,m. If 
the observed dust results from collisional grinding of larger solids, we infer mini- 
mum masses of the parent body population between 0.004 M e and 0.06 M®. We 
estimate grain production rates of ~10 10 g s _1 around A Boo and HR 1570; selec- 
tive accretion of this matter may help explain their peculiar surface abundances. 
There appear to be inner truncations in the dust clouds at 48 AU, 11 AU, 52 AU 
and 54 AU around HR 333, HR 506, HR 1082 and HR 3927, respectively. 
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1. INTRODUCTION 

The discovery with IRAS of a large infrared excess around Vega (Aumann et al. 1984) 
initiated the detailed study of other planetary systems. We can now use the dust emission 
from main sequence stars to constrain models for the origin and evolution of these environ- 
ments (Lagrange, Backman & Artymowicz 2000, Zuckerman 2001) with the long-term goal 
of developing a more comprehensive understanding of the formation and evolution of planets 
and related minor bodies. Most previous studies of the dust around main-sequence stars 
have been restricted to using broad-band infrared fluxes (e. g. Habing et al. 2001, Spangler 
et al. 2001). Here, we report measurements with the Infrared Spectrograph 13 (IRS, Houck 
et al. 2004) on the Spitzer Space Telescope (Werner et al. 2004) of 19 main-sequence stars. 

2. OBSERVATIONS 

We obtained spectra of main sequence stars with previously reported IRAS excesses with 
both the Short-Low (5.2 - 14.5 /im) and Long-Low (14.0 - 38.0 /im; A/AA ~ 90) modules 
on the IRS. Table 1 lists the six targets in our sample with strong infrared excesses. We did 
not detect a substantial excess around 13 other stars. An example of a star without a strong 
excess is HR 4732; its spectrum (from Spitzer aor 0003579392) is shown in Figure 1. The other 
stars without strong excesses are HD 16157 (aor 0003554560), HD 221354 (aor 0003565842), 
HR 818 (aor 0003554816), HR 1338 (aor 0003555584), HR 2015 (aor 0003556608), HR 3220 
(aor 0003557888), HR 3862 (aor 0003558400), HR 5447 (aor 0003560192), HR 8085 (aor 
0003564544), HR 8549 (aor 0003588608), HR 8799 (aor 0003565568) and Ross 128 (aor 
0003559168). 

In order to avoid time-consuming peak-up on our relatively bright targets with accu- 
rately known positions, we operated the observatory in IRS Spectral Mapping mode where 
a 2x3 raster (spatial x dispersion) centered on the star is performed. The slit positions were 
separated by half of the slit width in the dispersion direction, and by a third of the slit length 
in the spatial direction. 

We carried out the bulk of the reduction and analysis of our spectra with the IRS team's 
SMART program (Higdon et al 2004). We started with intermediate products from the 
Spitzer IRS data-reduction pipeline, that lacked only stray light and flat-fielding corrections. 
We extracted point-source spectra for each nod position using a uniformly weighted column 



13 The IRS was a collaborative venture between Cornell University and Ball Aerospace Corporation funded 
by NASA through the Jet Propulsion Laboratory and the Ames Research Center. 
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profile for which the variable width scales with the width of the instrumental pointspread 
function. We also made similar spectral extractions for observations of our photometric 
standard star, a Lac (A1V). Then we divided each target spectrum, nod position by nod 
position, by the spectrum of the standard star, multiplied the quotient by the appropriate 
template spectrum (Cohen et al 2003), and averaged the resulting spectra from corresponding 
nod positions. We plot the resulting spectra of our targets with IR excesses in Figure 1. 
Based upon comparisons of IRS spectra of non- variable calibration sources to ground-based, 
IRAS, and Spitzer IRAC flux densities, we estimate that the absolute spectrophotometric 
accuracy of these spectra is 10% (lo - )- It is important to distinguish this absolute photometric 
accuracy, which describes the vertical positions of the spectra on our plots relative to the 
flux density scale, from the point-to-point fluctuation in the spectra, which is a measure of 
the precision of our spectral feature strengths and therefore of our ability to identify spectral 
features and compare them to our models. In general, the accuracy of our spectral feature 
identifications and feature strengths is limited by the point-to-point scatter (currently much 
larger than the system noise, in bright objects) that is visible in the spectra. Therefore 
almost all features that significantly exceed the fluctuations, and are wider than the two- 
pixel IRS spectral resolution element, are significant. The lack of spectral features that we 
refer to in this paper is based on this criterion. 



3. INTERPRETATION AND ANALYSIS 

To interpret the data for the infrared excess, we need luminosities, L*, and ages, £*, of 
the stars. We take distances, D*, from Hipparcos, apparent magnitudes from the Yale Bright 
Star Catalog, assume no interstellar reddening, and use the bolometric corrections given by 
Flower (1996) to infer the total luminosities. For four stars, we adopt published estimated 
ages. For HR 506, the age is highly uncertain since Zuckerman & Song (2004) and Decin et 
al. (2000, 2003) estimate 0.3 Gyr and 3.0 Gyr, respectively. In view of the recently measured 
lithium abundance of 2.6 on the usual 12-point logarithmic scale (Israelian et al. 2004), we 
adopt the younger age. For HR 333 and HR 3927, we assume ages of 0.1 Gyr, representative 
of dusty main sequence A-type stars (Song et al. 2001). 



3.1. Grain Size 

The dust continua for A > 20 /zm show no spectral features. If the particles are spheres 
of diameter a, then any spectral features that might be present are weakened if a > \/tt 
(e. g. Spitzer 1978, Wolf & Hillenbrand 2003). Our data thus suggest that the diameter 
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of the grains is larger than 10 /im, as predicted by models such as those of Krivov, Mann 
& Krivova (2000) and found for the particles that produce the zodiacal light in the Solar 
System (Fixsen & Dwek 2002). 

3.2. Spatial Distribution and Dynamics 

A standard model of dust debris systems assumes that collisions among a population 
of large, unseen parent bodies produces smaller, detectable dust particles (Lagrange et al. 
2000, Zuckerman 2001). Some of the best studied systems such as HR 4796 A exhibit a 
well defined ring of material (Schneider et al. 1999). Very approximately, for circumstellar 
envelopes with L 30 /L* > 10~ 4 (where L 30 is the characteristic luminosity of the infrared 
excess in the 30 jum spectral window, or vL v at 30 /im), particle-particle collisions typically 
are the most important grain destruction mechanism (e. g. Krivov, Mann & Krivova 2000, 
Wyatt & Dent 2001). In these relatively dusty environments, the grains do not migrate 
far before they are fragmented into sufficiently small pieces that they are driven out of the 
system by radiation pressure. Thus, the spatial distribution of the dust reproduces the 
spatial distribution of the parent bodies which might be determined by the history of planet 
formation in the system. In contrast, in less dusty envelopes with L 30 /L* < 10~ 4 , mutual 
collisions are not so important, and the grains spiral inwards under the Poynting- Robertson 
drag. In such models, the spatial distribution of the grains need not reproduce that of the 
parent bodies. 

Motivated by these models of dust debris, we consider two possible fits to the reduced 
data. First, we fit the spectra by a function of the form: 

F v = K x v 2 + K 2 v~ x (1) 

where K x and K 2 are constants. The first term in equation (1) describes the photospheric 
emission while the second term represents the infrared excess. We show in Figure 1 the 
comparison of our best fits with the data, and we list in Table 1 our values for K 2 and 
L 30 /L* where L 30 is AnD 2 K 2 . Equation (1) fits the spectra for the two stars with the lowest 
values of L 30 /Z/*, HR 1570 and A Boo. However, the infrared excesses of four stars require 
another model since fits using equation (1) either underestimate the observed photospheric 
fluxes measured at A < 12 /zm (HR 333 and HR 3927) or the measured IRAS and/or ISO 
fluxes (Decin et al. 2000) at A = 60 /im (HR 506 and HR 1082). 

Our second method to model the infrared excess is to fit a single temperature to the 
dust particles so the total flux is fit with a function of the form: 

F v = C x v 2 + C 2 B u (T ex ) (2) 
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where B v is the Planck function. In principle, C\ must equal Ki, but the absolute level of 
our spectra is somewhat uncertain because of imprecise calibration. As can be seen in Figure 
1, we can fit our infrared spectra adequately with equation (2) with the values of T ex given 
in Table 1 ranging from 72 K to 114 K. 

The spectral energy distribution of the infrared excess can test the models. If Poynting- 
Robertson drag dominates, the spectrum varies as u^ 1 (e. g. Buitrago & Mediavilla 1985), 
the dust is smeared over a wide range of radii and equation (1) should fit the data. In 
contrast, in a model where collisions limit the lifetime of the particles, then inward drift is 
not as important and the dust is more concentrated near its birthplace. In environments 
with L 30 /L* > 10~ 4 , there might be a dust ring or at least a well defined inner hole to the 
dust distribution which achieves a characteristic maximum temperature. Within the limited 
IRS spectral range, equation (2) can fit the data. Inner dust truncations might be the result 
of the formation of low mass companions which might be planets (Kenyon & Bromley 2002), 
but other models can produce rings as well (Takeuchi & Artymowicz 2002). Also, since 
T ex is near 100 K, if water ice is an important constituent of the material, there may be 
substantial grain destruction by sublimation which contributes to the observed deficiency of 
warm grains (Jura et al. 1998). 

We now estimate some parameters for the circumstellar dust. We assume that the 
particles radiate like black bodies so that particles at temperature T lie at distance D from 
the star: 



For the two stars where the spectra can be fit by equation (1), there is no constraint on the 
maximum temperature and thus no evidence for any inner truncation in the dust distribution. 
There are four stars where we use equation (2) to fit the data, and we thus find a characteristic 
grain temperature, which, from equation (3), we can translate into a distance from the star. 
Using the values of T ex and given in Table 1, we find values for the inner boundary of 
the dust of 48 AU, 11 AU, 52 AU and 54 AU around HR 333, HR 506, HR 1082 and HR 
3927, respectively. Consistent with models for those dust debris systems where frequent 
collisions lead to smaller particles which are then ejected by radiation pressure, the stars 
with the strongest evidence for inner holes in the dust distribution have the highest values 
of L 30 /L*. If the inner holes are caused by inner planets, then in the future, with sufficiently 
high angular resolution, it may be possible to identify azimuthal asymmetries in the dust 
which might indicate gravitational perturbations by this unseen system of planets (e. g. 




(3) 



Wyatt 2003). 



We can also use our data to make a rough estimate of the minimum mass of the parent 
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bodies, Mp B - Since Poynting- Robertson drag provides a maximum grain lifetime, then in 
an approximate steady state, we can write (Chen & Jura 2001): 

M PB >^ (4) 

c 2 

This estimate of M PB is insensitive to the grain size and composition; the minimum values 
of Mp B derived from equation (4) are shown in Table 1. In systems with L^/L^ >> 10~ 4 , 
collisions among dust grains occur often, the Poynting-Robertson lifetime is thus much longer 
than the collision time, and we may severely underestimate M PB . There may also be large 
objects such as planets that are not experiencing rapid collisions and which do not contribute 
to the detected dust. 



4. DISCUSSION 

We show in Table 1 that the inferred lower bounds to M PB range from 0.004 M e to 
0.06 M e . The summed mass of the KBOs in the Solar System is uncertain and is currently 
estimated to be between 0.02 M e and 0.1 M© (Bernstein et al. 2003; Luu & Jewitt 2002). 
Typical analogs of the Kuiper Belt may have total masses < 0.1 M© (Jura 2004). Since we 
observed stars with notable infrared excesses, we probably selected objects with particularly 
massive systems of parent bodies. In the Solar System, the value of the infrared luminosity 
of the dust is ~ 10 -7 L (e. g. Fixsen & Dwek 2002). Either the mass of parent bodies 
around our sample of stars is much larger than the inferred minimum, or the collision rate 
of these objects is substantially greater than in the Kuiper Belt of the Solar System. 

Two of the stars in Table 1, A Boo itself and HR 1570, are A Boo stars, a subclass of 
~2% of all main-sequence A-type stars that is defined as having nearly solar abundances of 
C, N and O, but marked deficiencies of Fe elements by factors of 10-100 (Paunzen et al. 
2002). The explanation for this abundance pattern is not known. One theory to account for 
their surface abundances is that they have preferentially accreted some circumstellar matter 
(Venn & Lambert 1990), and within the star's atmosphere there is an element-dependent 
chemical separation controlled by gravitational settling and radiative acceleration at the base 
of the thin outer convection zone of the star (Turcotte & Charbonneau 1993). As shown in 
Jura (2004), the rate at which circumstellar dust is being produced scales as the infrared 
excess luminosity such that M ~ L 30 /c 2 where c is the speed of light. We therefore find 
from the results in Table 1 that M ~10 10 g s^ 1 for both A Boo and HR 1570. According to 
Turcotte & Charbonneau (1993), hydrogen-rich gas infall rates of >10 12 g s^ 1 with selective 
accretion of >10 10 g s~ x of metals in dust is required to explain the surface abundances of 
A Boo stars. While more sophisticated treatments are required to the accretion/diffusion 
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models, our results at least show that metal infall at the required rate appears to be occuring 
around A Boo stars. 

In models for the first 100 Myr of the Solar System, Petit, Morbidelli & Chambers (2001) 
have computed that approximately 50% of the ~5 M e of planetesimals initially within the 
asteroid belt had their orbits perturbed by Jupiter so they hit the Sun. Murray et al. 
(2001) have argued that such models can be generally extended to other main sequence 
stars. Assuming that this accretion of ~5 x 10 12 g s _1 of rocky material often occurs, 
the atmospheres of the A stars should be substantially polluted. We speculate that the 
main sequence stars with inner holes in their dust distributions may also be systems with 
asteroidal belts where gravitational perturbations by giant planets leads to a significant 
supply of metals to an atmosphere where they would be otherwise depleted. An infall of 
large rocks or asteroids may lead to observable variations in the abundance of the elements 
in the star's atmosphere (Cowley 1977) which could be different from the effects of dust 
inflows occuring around the A Boo stars. 

5. CONCLUSIONS 

We have obtained Spitzer mid-infrared spectra of 19 main-sequence stars. Our data are 
most useful for constraining the evolution of the solid material in regions between ~10 AU 
and ~60 AU from the stars. We conclude the following: 

• There are no strong spectral features; we suggest that the typical diameter of a dust 
particle is greater than 10 /xm. 

• The infrared excesses can be explained by dust arising from collisions between parent 
bodies whose minimum total mass lies between 0.004 M e and 0.06 M e . These derived 
minimum masses are comparable to the total mass of KBOs in the Solar System. 

• We estimate dust production rates of ~10 10 g s _1 around A Boo and HR 1570; selective 
accretion of this material may help explain their peculiar surface abundances. 

• There appear to be inner truncations in the dust clouds at 48 AU, 11 AU, 52 AU and 
54 AU around HR 333, HR 506, HR 1082 and HR 3927, respectively. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA through Contract Number 
1257184 issued by JPL/Caltech. F.K. is supported by the Spitzer Fellowship Program under 
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award Oil 808-001. Work at Cornell and the University of Rochester has been supported 
by JPL contract 960803 to Cornell and Cornell subcontract 31419-5714 to the University of 
Rochester. We are grateful to the entire team of scientists and engineers that have made 
Spitzer and the IRS function so well. 
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Table 1 - Stars with Excesses in the IRS data 



Star 


Sp. 


D 


L* 




K 2 




T 

- 1 ex 


M PB 


t(exp) 


AOR 3 






(pc) 


(L ) 


(Gyr) 




(io- 4 ) 


(K) 


(M e ) 


(s) 




HR 333 


A3V 


83 


34 


0.1 


0.16 


1 


97 


0.03 


24 


0003553536 


HR 506 


F9V 


17 


1.5 


0.3 1 


0.15 


0.9 


94 


0.004 


24 


0003553792 


HR 1082 


A3IV/V 


74 


12 


0.1 1 


0.15 


2 


72 


0.02 


24 


0003555328 


HR 1570 


AOV 


37 


15 


0.1 2 


0.14 


0.4 


88 


0.005 


24 


0003555840 


HR 3927 


AOV 


98 


44 


0.1 


0.21 


1 


98 


0.06 


24 


0003558656 


A Boo 


AOV 


30 


15 


0.2 2 


0.21 


0.4 


114 


0.01 


24 


0003559936 



1 Zuckerman & Song (2004); 2 Heiter, Weiss & Paunzen (2002). In this Table, K 2 (10 
Jy Hz) is defined in equation (1). 3 The Spitzer aor data file identification number. 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Fig. 1. — Plots of F v vs. A for the 6 stars in Table 1 plus HR 4732 which has, at most, 
a modest infrared excess. The fits from equation (1) are shown in red while the fits from 
equation (2) are shown in blue. We also display 3<r error bars at 25 /zm in green. The sum 
of the two terms in each equation is given by a solid line while dotted and dashed lines 
represent the photospheric fluxes and dust excesses, respectively. 



